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Andreev interferometry — the sensitivity of the tunneling current to spatial variations in the 
local superconducting order at an interface — is proposed as a probe of the spatial structure of the 
phase correlations in the pseudogap state of the cuprate superconductors. To demonstrate this idea 
theoretically, a simple tunneling model is considered, via which the tunneling current is related to 
the equilibrium phase-phase correlator in the pseudogap state. These considerations suggest that 
measurement of the low- voltage conductance through mesoscopic contacts of varying areas provides 
a scheme for accessing phase-phase correlation information. For illustrative purposes, quantitative 
predictions are made for a model of the pseudogap state in which the phase (but not the amplitude) 
of the superconducting order varies randomly, and does so with correlations consistent with certain 
proposed pictures of the pseudogap state. 
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I. INTRODUCTION 



A range of experimental investigations have indi- 
cated that underdoped high-temperature superconduc- 
tors (HTSCs) exhibit intriguing properties at tempera- 
tures above the superconducting transition temperature 
T c . Most notably, these materials show a strong sup- 
pression in the single-particle electronic spectral weight 
at low energies, even at temperatures far above T c 
a property referred to as the pseudogap. A number of 
scenarios have been proposed to account for this loss of 
spectral weight p|-pT|, several of which invoke the notion 
that remnants of superconductin g co rrelations remain in 
the non-superconducting state i.e., that pairing is 

established locally but that it lacks the long-range coher- 
ence in phase necessary for true superconductivity. 

To make progress with understanding the nature of 
the pseudogap regime, having experimental access to the 
spatial structure of the correlated electronic state is likely 
to be of considerable value |l2| . The aim of the present 
Paper is to identify one possible scheme, involving low- 
voltage mesoscopic conductance measurements, for prob- 
ing this structure experimentally, and to describe this 
scheme within the context of a simple theoretical model. 

The basic idea is this. Let us adopt as a work- 
ing hypothesis the picture of the pseudogap regime in 
which superconductivity is established locally, but in 
which the presence and motion of vortices in the su- 
perconducting order parameter cause the phase of the 
superconducting order parameter to be randomized be- 
yond certain correlation length- and time-scales Jl5| . 
The effects of such phase fluctuations on the single- 
particle properties of underdoped cuprates have been ex- 
plored in Rcfs. Jl^Jl^]. Now, the low- voltage conduc- 
tance of a normal-to-superconducting junction includes 
contributions associated with the Andreev reflection of 



quasiparticles from the superconducting condensate [ fl8[ . 
What about the low-voltage conductance of a normal- 
to-pseudogap junction? Given the picture of the pseudo- 
gap regime outlined above, and assuming that tunneling 
through the junction occurs on a time-scale faster than 
the time-scale for vortex rearrangement, we anticipate 
that there will be contributions to the conductance due to 
the Andreev reflection of quasiparticles from the local su- 
perconductivity. However, owing to the phase sensitivity 
of the Andreev reflection process |^] , any spatial varia- 
tion in the phase of the superconducting order parameter 
over the junction would tend to cause (diffraction-like) in- 
terference of the quasiparticle/hole waves that have been 
Andreev reflected from the junction, and thereby dimin- 
ish the associated contribution to the conductance. 

Now suppose that the normal contact in a normal-to- 
pseudogap junction has a characteristic linear dimension 
L. If L is smaller than the characteristic phase-phase 
correlation (i.e. inter-vortex) length £^ (e.g. the smaller 
contact in Fig. Q) then, at any instant, rather little phase 
variation would be expected over the contact and the An- 
dreev contribution to the conductance should be barely 
diminished. However, if L is substantially larger than £^ 
(e.g. the larger contact in Fig. |l|) then considerable phase 
variation is expected over the contact, and the Andreev 
contribution to the conductance is likely to be strongly 
suppressed. Measurements made using a range of meso- 
scopic contact sizes thus have the capability of providing 
a direct probe of the spatial correlations of the phase of 
the superconducting order parameter at various temper- 
atures within the pseudogap regime. 

Let us emphasize that the concept of Andreev interfer- 
ometry is by no means new; indeed, several groups have 
considered this concept both theoretically |2!| and ex- 
perimentally p4|] . However, to the best of our knowledge 
this interferometry has primarily been considered in con- 
texts in which reflection is from a truly superconducting 
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region (rather than from a pseudogap region) , and in set- 
tings in which the phase has an average value that varies 
in a relatively simple way in space (such as on either side 
of a Josephson junction). Here, we are considering a set- 
ting in which the interferometry is being used as a probe 
of the superconducting fluctuations. 




FIG. 1. Schematic depiction of an instantaneous configu- 
ration of the superconducting phase in the pseudogap state 
(arrows). The two shades of gray indicate two possible (nor- 
mal-state) contact areas on the pseudogap (i.e. white) sub- 
strate. Whereas the smaller (i.e. darker) contact abuts a re- 
gion of nearly uniform phase, the larger (i.e. lighter) contact 
exhibits regions of considerably differing phase. 

It should be mentioned that in recent work Choi et 
al. |2j| considered the issue of whether or not the zero- 
bias tunneling conductance peak pq ] would survive at 
temperatures above T c . This work involves applying the 
BTK technique to the physical picture of the pseu- 
dogap regime explored, e.g., in Refs. |^6|,0. It amounts 
to a computation of the conductance of a normal-to- 
superconductor interface for a d-wave superconductor in 
a uniform supercurrent-carrying state, this conductance 
then being averaged over a Gaussian distribution of uni- 
form supercurrents (in order to model the pseudogap 
state). This yields a conductance dependent upon the 
statistical distribution of local values of the supercur- 
rent arising from varying vortex locations. In contrast, 
the present work focuses on the spatial correlations of 
the phase in the pseudogap regime and, specifically, how 
such correlations may be accessed experimentally. 



the tunneling formalism, and show how this conductance 
depends on the pseudogap phase-phase correlation func- 
tion. To this end, we adopt as the tunneling Hamiltonian 
i?T @1 between a normal state (N) and a pseudogap 
state (P): 
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where the position s lies on the normal side of the junc- 
tion and the position r lies on the pseudogap side. The 
operators c s a (or c\ a ) and d r a (or d\ a ) respectively an- 
nihilate (or create) quasiparticles with spin projection a 
on the normal side at s and on the pseudogap side at r. 
We choose the interface to be in the plane z = (where 
{x,y,z} are Cartesian coordinates and { e z } are 

the corresponding basis vectors) and, accordingly, de- 
compose vectors such as s and r into components parallel 
(e.g. <r) and perpendicular (e.g. s z e z ) to the interface so 
that s = cr + s z e z and r = p+r z e z . This choice, together 
with the assumption that tunneling only occurs locally at 
the interface leads us to assert that the tunneling matrix 
elements i rjS are given by 



t r , s =toaS^(s z )5^(r z )S^(a-p), 



(2.2) 



where a is a microscopic length scale characterizing the 
thickness of the "active" layer for tunneling of particles 
and to is the typical energy scale for this process. 

Next, we compute the current /(V) as a function of 
the voltage V. To do this, we consider the expectation 
value of the tunneling current operator [— eQ^, Ht]/i?i, 
where — e is the electron charge: 
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with respect to the full Hamiltonian for the system, i.e., 
H = _Hn + Hp + Ht, where H^/p is the Hamiltonian 
for the normal/pseudogap side and — cQn is the charge 
operator for the normal side. In fact, it is convenient 
to obtain I(V) perturbatively in the tunneling ampli- 
tude toa by applying the Matsubara technique to the 
imaginary-time dependent tunneling current I{t{) |p8| . 
The lowest-order term, which is of order |£oo| 2 , represents 
the normal (i.e. single-particle) current. This contribu- 
tion is suppressed at low voltages due to the presence of 
a gap at low energies on the pseudogap side. The next- 
order contribution to /, which is of fourth order in toa, 
is given by 



II. TUNNELING CURRENT FOR A 
NORMAL-TO-PSEUDOGAP JUNCTION 

We now illustrate the ideas of Sec. | by computing the 
conductance of a normal-to-pseudogap junction within 
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where (• ■ -)p/n indicates an equilibrium expectation value 
with respect to ifp/N, and measures the inverse 
temperature. Operators such as c ro .(r) [or d a (t)] 
are interaction-picture operators, i.e., e K ° T c r a cT KaT (or 
e K ° T d s IT e~ KoT ), where K = iJ N - ^n<2n (or K = 
Hp — j-ipQp). Here, /in (or fip) is the chemical potential 
on the normal (or pseudogap) side, and Qn (or Qp) is 
the charge operator for the normal (or pseudogap) side. 
The physical current I(V) is given by the imaginary part 
of I(t±) after making the following analytical continua- 
tions: iuig — > i0 + , «H — > eV (i.e. the voltage across the 
junction), and — ih/rx — > t (i.e. the time). 



To apply Eq. (2.4) to the setting at hand, namely one 
side of the junction being normal and the other being 
in the pseudogap regime, we shall need to evaluate the 
two two-particle Green function factors that feature in 
it, one for the normal side and one for the pseudogap 
side. For the normal side, we assume that the corre- 
sponding two-particle Green function is factorizable into 
two single-particle Green functions, i.e., 

= G N ( Pl , p 4 ; n - r 4 ) G N (p 2 , p 3 ; r 2 - r 3 ), (2.5) 

where G N (p, p'; r - r') = (T T c p Jr) c\, Jt')) is the 
single-particle Green function on the normal side. On 
the other hand, for the pseudogap side we adopt a model 
in which the pseudogap state is a superconductor that 
is "disordered" by a static pattern of vortices and char- 
acterized by a static phase-phase correlator. We do not 
require information concerning the dynamic phase-phase 
correlation function because we are assuming that the 
Andreev reflection process is rapid, compared with the 
time needed for vortices to substantially rearrange the 
phase structure. To support this assumption, let us note 
that the time-scale associated with Andreev reflection 
is of order tar ~ £,h/vf, where vp is the Fermi veloc- 
ity of the incoming electron and £ a is the amplitude- 
fluctuation correlation length (i.e. the Cooper-pair size) 
on the pseudogap side. Then, by using the estimates 
vp « 10 7 cm/s and £ a « 1 nm (typical for a HTSC) we 
find that tar ~ 10~ 14 s. The experiments of Corson et 
al. [ p9[ indicate that the vortex-pattern rearrangement 
time corresponds to frequencies in the terahertz range 
(i.e. is of order 10 _9 s) so that, at least as a starting 
point, we may neglect the dynamics of the vortices. Thus, 
we assume that quasiparticles incident from the normal 
side effectively encounter, and are Andreev reflected by, 
a static pair-potential that has a nonzero amplitude (ex- 
cept at the vortex cores, which are small) and a spatially 
random phase. With this in mind, we characterize the 
pseudogap side by the anomalous Green function 

F p (r,r';r,r')= <T T d ra (T)d r ,, T (r'))p (2.6) 
= / (r-r';r-T')e^W + ^ r '»/ 2 , 

where the phase </>(r) varies randomly in space, and the 



function /o is given by the value it takes in a conventional 
superconductor Q, i.e., 

/o(r;r)= /"tjSU £ e-V k ^ 



71 — — OC 

A k \ ( 1 



1 



2Eu J \ iuj n + Eu iu n - El 



(2.7) 



Here, E k [= + A£, with £ k = (h 2 k 2 /2m) - n] is the 
excitation energy in the pseudogap material and is 
the gap amplitude. The Matsubara frequencies cu n are 
defined to be w„ = (2n + l)ir/fl for integer n. For the 
sake of simplicity, we now focus on the case of s-wave 
pairing and thus set Ak = A |5l| . We approxim ate the 
two-particle Green function on the P side in Eq. (2.4) by 
making a Gorkov factorization into the anomalous Green 
function F p and its conjugate. In principle, there will 
also be a contribution associated with factorization into 
normal Green functions. However, these contributions 
are suppressed at low voltages, and thus it is adequate 
to take for the current 
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X /o(Pl - P2'i iuj n) /o(P3 - Pi> iuJ n + 2iuJ S ), (2-8) 

where the limit A indicates that the interface integrals 
are constrained to the area of the contact between the 



N and P regions. Equation (2.8) may be expressed dia- 
grammatically, as shown in Fig. 2, where one-arrow lines 
denote normal Green functions, two-arrow lines denote 
anomalous Green functions (and the straight line denotes 
the interface). We see that this contribution to the cur- 
rent involves the correlation of an electron and a hole 
propagating on the normal side, mediated by the static 
random pair-potential on the pseudogap side of the junc- 
tion. 




FIG. 2. Diagram depicting the leading-order (in tunneling 
amplitude) phase-sensitive contribution to the current. 
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Equation which represents the leading contribu- 

tion due to Andreev reflection at an interface, may be 
considerably simplified in situations in which £ a <C £0 
(i.e. the phase order which we are interested in probing 
via Andreev reflection persists over length scales much 
larger than the pair size), which is not only the case for 
the usual NS interface setting but also for the present 
NP setting. To support this assertion, we obtain the ra- 
tio of these two length-scales by examining the results of 
the Corson et al. experiments on the high-frequency AC 
conductivity of Bi2Sr2CaCu2 08+5 p9[ ). According to the 
analysis of Corson et al., in the pseudogap state the ratio 
of £0 to £ a is related to the vortex diffusion time t via 
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(2.9) 



where Oo is a parameter determined by Corson et al. to 
be 1.14 x 10 14 s _1 . From Fig. (4) of Corson et al., we see 



that at T = 75 K, r ~ 1CT 12 s, 



so that ^/£a 



20. 



The significance of the separation of the length-scales 
£ a and ^ in the present context follows from the fact tha t 
the function /o has spatial range £ a - Thus, in Eq. (2.?) 



the spatial integrations over the coordinates {Pi}i—% may 
be simplified because /o varies rather more rapidly in 
space than do the other factors in the integrand. This 
allows us to make the approximation 
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!%/ (fc z e z ,zc), (2.10) 

where fo is the (three-dimensional) Fourier transform of 
/o and k z is the momentum component perpendicular to 
the interface. By using this approximation, we obtain 
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Having derived Eq. (2.11), an equation applicable to any 



given realization of the phase field 4>(p) on the P-side of 
the interface, we conclude the present section by perform- 
ing the averaging this current over an as-yet-unspecified 
distribution of phase fields. As discussed above, the time- 
scale for the tunneling process is shorter than the time- 
scale for phase rearrangement. Thus, it is appropriate 
to proceed as we have, by first computing the current 
for a fixed realization of the phase field, and then to 
construct the time- averaged current, averaged over times 
longer than the phase rearrangement time, by averaging 
the current over an appropriate (in this case, equilibrium) 
distribution of phase fields. Denoting such averaging by 



] , and introducing the appropriate phase-phase corre- 



lator 



(2.12) 



we arrive at a for mula for the time-averaged current 
[/(tjl)], i.e., Eq. ( [2.1l| ) but with the phase factors 
expi(^(p x ) - (j)(p 2 )) replaced by g(p 1 - p 2 ). For con- 
venience, we express the normal-side Green function in 
terms the corresponding spectral function A: 



GN (Pi ; P2;* w r, 



de A{p 1 ,p 2 ;e) 
27T iuj n — e 



(2.13) 



Then we may perform the integrations over k z and k z (by 
converting them to energy integrals) , as well as the sum- 
mation over Matsubara frequencies. By performing the 
necessary analytic continuations and taking the imagi- 
nary part, we obtain an expression for the tunneling cur- 
rent I(V) through a mesoscopic interface between a nor- 
mal metal and a material in the pseudogap state: 



l(V) = Jj 2 P iJji 2 P2 g(p 1 -p 2 ) 

:{n{e-2eV)-n{e)} 
(2.14) 
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where n(e) = (exp(/3e) + 1) _1 is the Fermi distribu- 
tion function and i>p [= m/ (27rh 2 kp)~ 1 with kp being 
the Fermi wave- vector] is the one-dimensional density of 
states on the P-side. 



III. CASE OF CLEAN NORMAL-METAL 
CONTACT 

A. General considerations 



In this section we pursue the evaluation of Eq. (2.14) 
for the case of a normal contact that is perfectly clean, 
in the sense that the spectral function A c (with the su- 
perscript C standing for clean) has the form appropriate 
for a pure metal: 



A c (p;e) = 27rS^(e p -e), 



(3.1) 



in which e p = p 2 /2m — p. The (three-dimensional) 
Fourier transform of this quantity is given by 

A c (r,r'; £ )^ J ^P A c (p;e)e i ^- r 'y\ (3.2a) 



{2iTh) s 

m sin{^2mh' 2 {e + p)\r- r'|} 
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By inserting this expression into Eq. (2.14), and limit- 
ing our attention to low temperatures (i.e. k-^T <C eV, 
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with kft being Boltzmann's constant) and low voltages 
(i.e. eV <C A) we obtain an equation for the low- 
voltage conductance as a functional of the pseudogap 
phase-phase correlation function 
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B. Illustrative example: BKT correlations 



The main conclusion of Sec. Ill A is that the contri 



bution of Andreev reflection to the tunneling current is 
sensitive to spatial inhomogeneity of the superconduct- 
ing phase, such as is proposed to exist in the pseudogap 
state. For the purposes of illustration, we now exam- 
ine a specific example of how the current enhancement 
due to Andreev reflection is increasingly suppressed, with 
increasing area, due to destructive interference. In this 
example, we assume that the phase-phase correlations in 
the pseudogap state are adequately modeled by those as- 
sociated with the Berezinskii-Kosterlitz-Thouless (BKT) 
theory of the two-dimensional XY model [^3|- ^5|. The 
relevance of this theory to the cuprate materials ]36|] orig- 
inates in the fact that their pronounced planar charac- 
ter causes the intermediate length-scale electronic struc- 
ture to be characterized by two-dimensional XY behav- 
ior, which is expected to cross over to three-dimensional 
XY behavior only very close to the transition. In order 
to compute the current for this BKT scenario, we need 
a form for g(p). On length-scales short compared with 
the phase-phase correlation length the function g{p) 
approaches unity; on length-scales long compared with 
£0, it exhibits exponential behavior l3j|. As we are only 
seeking an illustrative computation of the current, the 
exact details of this crossover are unimportant, and thus 
we adopt the form 
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and we take the interface to have the shape of a disk 



of radius L. Inserting Eq. (3.4) into Eq. (3.3), we see 
that the low-voltage Andreev conductance per unit area 
through the interface has the form 
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(3.5c) 



Here, the subscript 1 indicates that the integrals are 
taken over disks of unit radius. The prefactor is the 
limiting value of the conductance per unit area in the 
limit of large interface area. [Note that Tq vanishes for 
the case of no phase coherence (i.e. kp^ = 0).] 

One of our primary concerns is how the varying of 
the interface size would provide information regarding 
the structure of the phase correlations; this information 
is contained in the function fc, which depends only on 
the dimensionless quantities kp^ and L/£$. For generic 
values of its arguments, the form of f c can be deter- 
mined only via numerical integration; however, its be- 
havior can be determined in various physically relevant 
asymptotic limits. To begin with, let us assume that the 
phase correlations persist over length-scales that are long 
compared with the Fermi wavelength on the normal side 
(i.e. /cf£<^> ^> 1), and let us consider varying the interface 
size. For small interface sizes (i.e. fcp 1 <C L <C ^), fc in- 
creases logarithmically with L (i.e. fc ~ In k-pL/ In kp^). 
In the opposite regime of large interface sizes, we expect 
that Andreev reflection will occur from independent "do- 
mains" of uniform phase (so that, e.g., the doubling of 
the area should double the conductance). Indeed, for 
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for any value of kp^. 
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Dimensionless contact size, L/t,^ 

FIG. 3. 1 — fc (i.e. the departure of the dimensionless An- 
dreev conductance per unit area from its large-area limit) as 
a function of the dimensionless size of the interface L/^, for 
the case of = 100 (computed numerically). For L much 
smaller than the zero- voltage conductance per unit area 
is much smaller than its asymptotic value. 

To study the behavior of fc for intermediate values of 
L, we perform the integrals in Eq. ( 3.5c| ) numerically. In 
Fig. H| we show 1— fc as a function of L/£$ for the case of 
long-range phase correlations (i.e. kp^ — 100). Thus, as 
discussed Sec. |, a series of mesoscopic conductance mea- 
surements involving a range of contact sizes is expected 
to be rather sensitive to the characteristic length-scale of 
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phase correlations in the pseudogap state. 

IV. CASE OF DISORDERED NORMAL-METAL 
CONTACT 

A. General considerations 



In Sec. [II we investigated the conductance of a meso- 



scopic normal metal-to-pseudogap junction for the case 
of a perfectly clean normal metal. We now address the is- 
sue of the sensitivity of the main result (i.e. that this con- 
ductance contains information regarding the spatial ex- 
tent of the pseudogap-side phase correlations) to the as- 
sumption that the contact is a per fect ly clean metal |37| . 
Specifically, we examine how Eq. (3.3) is modified by the 
presence of disorder in the normal-metal contact. As we 
shall see, in the presence of disorder the most significant 
contribution to the conductance is associated with so- 
called Cooperon diagrams, familiar from the theory of 
the weak-localization corrections to the conductivity of a 
disordered metal fl38 |. 

As is conventional, we take the disorder to be due to 
uncorrelated point-like impurities, which scatter the elec- 
trons elastically. Moreover, we assume that the dephas- 
ing length is long, compared to both £^ and the mean 
free path I (which characterizes the strength of the disor- 
der and is related to the scattering time r via t = vft). 
Although we are focusing on situations in which is 
larger than the interface size so that one expects 
substantial sample-to-samplc fluctuations (which may in 
fact be interesting to study), we shall restrict our atten- 
tion solely to the disorder av erage of the current. Then, 
averaging the current in Eq. ( 2.14 ) over configurations of 
the potential scatterers on the normal side, an averaging 
that we indicate via (■ ■ -)dis, we arrive at 

(J(VO)dto = ~M 4 *3 J d 2 Pl J d 2 p 2 g{ Pl -p 2 ) 

r 2eV+)i £2 



& A 2_ (eF _ e)2 H £ - 2e ^)-^)} 

x (A D ( Pl , p 2 ; e) A D ( Pl ,p 2 ; 2eV - e)) dis , (4.1) 

where the superscript D refers to the disordered case. 
The disorder-averaged product of spectral functions 
(A D (p 1 , p 2 \ e) A B (p 1 , p 2 ; 2eV — e))dis contains contribu- 
tions that extend over length-scales | p x — /o 2 1 much larger 
than I. These Cooperon contributions provide the mech- 
anism for the transmission of the phase-sensitive informa- 
tion that would be probed in Andreev interferometry ex- 
periments involving a disordered normal- metal contact. 



B. Semiclassical picture 

Before deriving our result for the contribution of the 
Cooperon to (I(V))di S , we pause to motivate physically 



why this particular contribution is significant. In the con- 
text of weak localization, the Cooperon contribution to 
the conductance is usefully pictured in terms of construc- 
tive interference of pairs of paths involving the scattering 
of electrons from impurities in reverse order Efl] . This 
interference tends to "localize" electrons, thus causing a 
reduction in conductivity. In the present context, how- 
ever, the origin of the Cooperon is slightly different. To 
see this, consider the amplitude A%£ for an electron at 
position x in the normal metal to scatter from a sequence 
of impurities labeled by the index a, to Andreev reflect 
at the position r on the interface, to then scatter from 
the sequence of impurities labeled (3, and finally to return 
to the position x. Then, the probability for an electron 
leaving x and reflecting from the interface to return to x 
as a hole is given by the squared modulus of the sum of 
such amplitudes, i.e., 



0'a' 



(4.2) 
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FIG. 4. Schematic depiction of two semiclassical trajecto- 
ries in which an electron leaves position x in the normal region 
(N), undergoes multiple elastic scattering events, then under- 
goes Andreev reflection at the NS interface (horizontal line), 
and then returns to x as a hole via the same scatterers but 
in reverse order. Full (dashed) lines represent electron (hole) 
trajectories; crosses represent impurity scattering potentials. 

As is well known, the amplitude depends sensi- 
tively on the specific real izat ion of the disorder; thus, 
the right-hand side of Eq. (4.2) contains many terms that 
are disorder-dependent complex numbers. These contri- 
butions to -P(x) average to zero upon disorder averaging. 
However, amongst the collection of amplitudes there is a 
special subset describing processes in which the hole, as 
it returns from the interface to x, does so via the same 
collection of impurities visited by the electron on the out- 
bound segment of the trajectory but in the reverse order. 
If we denote the reverse of the sequence of impurities (3 



G 



by the sequence f3 then this special subset consists of the 
amplitudes ; these amplitudes have the form of real 
numbers, regardless of the specific locations of the impu- 
rities, except for a factor due to the phase shift associated 
with Andreev reflection. To see this, consider, e.g., the 
left-hand pair of paths (electron and hole) in Fig. [|. The 
electron path (full line) originates at x, scatters from im- 
purities at positions 1, 2 and 3, and then Andreev reflects 
as a hole. The hole then propagates back to x, scattering 
from the impurities at positions 3, 2 and 1 before return- 
ing to the position x. The dynamical phase acquired by 
the electron as it propagates to the interface is canceled 
by a phase of the opposite sign acquired by the hole. 
Thus, the amplitude for an electron at x to return as a 
hole at x depends only on the phase of the condensate at 



\ aa 
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(4.3) 



Thus , on e sees that the most significant contribution to 
Eq. (4.2) is given approximately by 



To analyze (I(V))dis we make two further simplifying 
assumptions. First, as in the clean case, we limit our at- 
tention to low temperatures (i.e. k^T <C eV), and thus 
we obtain 
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(4.7) 



where we have inserted the explicit real-space expression 
for the Cooperon. Second, by making the restrictio n to 
low voltages (i.e. eV <C A) |52), we may in Eq. (4.7) 
replace A 2 /(A 2 — (eV — e) 2 ) by 1. Furthermore, in 
the presence of disorder one has the natural length-scale 
Ly = WhD/eV. At sufficiently low voltages and small 
interface sizes, Ly will be much larger than typical val- 
ues of I Pi — p 2 |j so that one may expand to lowest order 
in L/Ly, thus obtaining 



P(x) 



V^ e i(<|i(r)-^(r')) ) 



(4.4) (I(V)) dh 



and hence that -P(x) is sensitive to the nature of the 
pseudogap phase-phase correlations, a sensitivity similar 
to that embodied in Eq. (3.3). 
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C. Microscopic calculation 

The explicit computation of the contribution of the 
Cooperon directly follows the usual analysis found in the 
context of weak localization; following Rammer |38|, we 
find that the disorder-averaged product of spectral func- 
tions has the form 



(A D (x,x', £ )A D (x,x', £ ')) di 



(lis 



(4.5) 



2tt 
n 



d 3 Q 



. e iQ.(*-xO/» {(7(Q ; e _ e ')+C(Q, e' - e)} , 



where the Cooperon propagator C(Q,w) = {—iuu/h + 
DCph^ 2 )^ 1 , the diffusion constant D = w|r/3 (in three 
dimensions), and ^n[= k-pm/(2iT 2 h 2 )] is the n orm al-side 
density of states. Inserting Eq. (fT^) into Eq. (LI) leads 
to the expression 
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x{C( Pl - p 2 , 2(eV - e)) + C( Pl - p 2 , 2(e - eV)) } , (4.6) 

where C (p, e) is the (three-dimensional) Fourier trans- 
form of C(Q,e). 



where Lp = ^JfiD /cf- Thus, as found in Sec. Ill for 
the case of a clean normal metal contact, the low- voltage 
conductance of a disordered metal-to-pseudogap junc- 
tion also contains information regarding the pseudogap 
phase-phase correlation function. 



D. Illustrative example: BKT correlations 



In this section we examine the area-dependence of the 
low-temperature and low-voltage conductance of a disor- 
dered normal metal-to-pseudogap junction for the case 
of BKT correlations. As in Sec. H, we assume that 
phase correlations decay in an exponential fashion, con- 
sistc nt w ith the BKT scenario. Our starting point is thus 
Eq. (4.8), tog ethe r with the model of phase correlations 
given b y Eq . (3.4). By considering the V — > + behavior 
of Eq. (4.8) we arrive at the low-temperature conduc- 
tance per unit area (for the case of an interface having 
the shape of a disk of radius L): 
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(4.9b) 
(4.9c) 



7 



1.0 
^ 0.8 
0.6 
0.4 
0.2 



2 4 6 8 10 

Dimensionless contact size, L/^ 

FIG. 5. 1 — /d (i.e. the departure of the dimensionless An- 
dreev conductance per unit area from its large-area limit) as 
a function of the dimensionless size of the interface Lj^ for 
the case of a disordered normal-metal contact (computed nu- 
merically). For L much smaller than the zero- voltage con- 
ductance per unit area is much smaller than its asymptotic 
value. 



By evaluating the integrals in Eq. ( 4.9c ), we obtain 



/d = 1 



Li 
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(4.10) 



where I\ is a modified Bessel function and L\ is a mod- 
ified Struve function. The asymptotic behavior of /d is 
linear for small L (i.e. /d ~ (8/37r)(L/^) for L/^ <C 1); 
for large L it approaches unity as an inverse power of L 
[i.e. / D « 1 - (2^/ttZ) for Lfa > 1]. In Fig. | we show 
how this function crosses over between these two limits. 
As with the case of the clean contact, the conductance 
shows marked sensitivity to the phase-phase correlations 
of the pseudogap state. 



V. CONCLUDING REMARKS 

In this Paper, we have proposed and explored theo- 
retically the possibility of using Andreev interferometry 
to probe the spatial structure of the phase correlations 
in the pseudogap state of the cuprate superconductors. 
The viability of this technique rests on the sensitivity 
of the tunneling current across mesoscopic normal-to- 
pseudogap junctions to spatial variations in the local su- 
perconducting order in the pseudogap state, as well as 
the possibility of using junctions having a range of areas. 

By considering a simple tunneling model, we have es- 
tablished a relationship between the tunneling current 
and the equilibrium phase-phase correlator characteriz- 
ing the pseudogap state. We have considered the cases in 
which the normal region (i.e. the contact) is either a clean 
or a disordered metal. In both cases, we have assumed 
that phase coherence length for quasiparticles on the nor- 
mal side is greater than the contact size. If this condition 
is not met then, is not the case then, throughout our re- 



sults, the interface size would have to be replaced by the 
dephasing length. 

To illustrate this Andreev interferometry proposal, we 
have applied our general results to a simple model of the 
pseudogap phase-phase correlations, which is intended to 
mimic the BKT correlations relevant to certain proposed 
pictures of the pseudogap state. Our considerations sug- 
gest that measurements of the low-voltage conductance 
of mesoscopic tunnel junctions of varying areas between 
normal-state and pseudogap-state regions would reveal 
information about the phase-phase correlations in the 
pseudogap state. 
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